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We show that radio-frequency spectroscopy on weakly-bound molecules is a powerful and sensi- 
tive tool to probe molecular energy structure as well as atomic scattering properties. An analytic 
expression of the rf excitation lineshape is derived, which in general contains a bound-free compo- 
nent and a bound-bound component. In particular, we show that the bound-free process strongly 
depends on the sign of the scattering length in the outgoing channel and acquires a Fano-type profile 
near a Feshbach resonance. The derived lineshapes provide an excellent fit to both the numerical 
calculation and the experimental measurements. 
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I. INTRODUCTION 

Radio-frequency (RF) spectroscopy is widely applied 
to many experiments on atoms or molecules, for which 
an exquisite energy resolution can be achieved. This is 
due to the magnetic coupling nature of the RF transition 
and the associated long coherence time of the systems. In 
recent experiments on ultracold weakly-bound molecules, 
rf spectroscopy also allows a precise determination of the 
molecular binding energy P| and the pairing gap in a de- 
generate Fermi gas 0. In the latter case, the high energy 
resolution reveals the fermionic nature of the pairing in 
the Bardeen-Cooper-Schrieffer (BCS) regime. 

In this article we investigate theoretically the radio- 
frequency excitation spectrum of ultracold molecules. 
Previous experimental work in this regime shows that 
molecules are dissociated upon receiving the RF photons 
P, 0. The associated bound-free excitation lineshape 
is characteristically broader and highly asymmetric com- 
pared to that for atoms. Here, we provide an intuitive 
picture to model and derive a simple analytic formula for 
the excitation rates and the lineshapes. We show that 
much information regarding cold collision properties can 
be extracted from the lineshape. In particular, the spec- 
trum is very sensitive to tuning a Feshbach resonance in 
the final state scattering continuum, where the rf spectra 
will evolve into two components: a bound-bound transi- 
tion to the newly formed molecular state and a remaining 
weaker bound-free transition. 

Our analytic results provide excellent fits to both the 
numerical calculations and the experimental measure- 
ments from the Innsbruck group 0, Q . We argue that 
radio-frequency spectroscopy on cold molecules has sev- 
eral advantages in probing ultracold collision properties 
over the conventional collision measurements. These in- 
clude a high energy resolution, the flexibility to probe dif- 
ferent channels and its insensitivity to the sample density 
and temperature. 

In this paper, we first introduce our model and derive 
the lineshape (Sec. II). We then compare our results to 
the numerical calculation (Sec. Ill) and discuss the spec- 



tral feature (Sec. IV). Finally, we compare our results to 
the experimental data (Sec. IV). 



II. MODEL 



We consider a weakly-bound molecule in the state \m). 
The bound state consists of two atoms with binding en- 
ergy Eb = tl? /2a{a — rg)^ > relative to the dissociation 
continuum A Here h is Planck's constant, ^ is the re- 
duced mass of the two atoms, a is the scattering length in 
the scattering channel A and rg is the interaction range 
of the van der Waal potential, which varies with inter- 
atomic separation r as — Cg/r^: 
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For weakly-bound molecules, we assume a 3> rg. The 
molecule is initially at rest and a radio-frequency pho- 
ton with energy Enp couples the molecule to a different 
channel A', characterized by the scattering length a'. For 
a' < 0, the final state is a continuum and the excited 
molecule dissociates; for a' > 0, a stable bound state 
|m') is also available in the final state channel and the 
molecule can either dissociate or be driven to the bound 
state |m'). In the following, we assume that both scat- 
tering channels are in the threshold regime (a, \a'\ ^ rg). 
Atoms, bound or unbound, then have the same internal 
wavefunction. We also assume the relative energy of the 
continuum threshold to be Ea' — Ea = Eq > 0. 



A. Bound-free transition 

First, we consider the bound- free transition. Energy 
conservation gives 



Erf =Eo + Eb + K, 



(2) 



where K — h^k^ /2ii > is the kinetic energy of the 
outgoing wave and k is the associated wavenumber. For 
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EjjF < Eq + Eh, the transition is forbidden. For Ej^p > 
Eq + Eh, the bound-free transition rate from the initial 
state |m) to the final state \K) is given by the Fermi's 
golden rule: 



extreme situations including a = a' , where Ef vanishes 
and a' = ±00, where the Wigner threshold law fails, will 
be discussed in later sections. 

The integrated rf line strength Sf is given by 
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where 



FfiK) 



i^*K{r)(l>m{r)dr 
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is the bound-free Franck-Condon factor per unit energy, 
Ml/ 2 is the RF interaction for the rf Rabi frequency fJ, 
4>m{T) is the bound molecular wavefunction in channel 
A, and (f) is the energy-normalized s-wave scattering 
wavefunction in channel A' . 

We can evaluate the bound-free Franck-Condon factor 
Ff based on the asymptotic behavior of the wavefunc- 
tions. 
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where 8' is the scattering phase shift in channel A'. Com- 
bining Eqs. lO and ©, we have 
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Since k <C l/j'Oj we can use the low energy expansion 
of the scattering phase shift , 



kcotS' = + ^fc2 + 0(/fc'') , 



(8) 



where the effective range depends on the scattering 
length a' and the interaction range tq. The explicit form 
of Tg for a van der Waals potential is given in Ref. in 



the limit 



a'\ > ro: 
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By taking the leading term in the expansion in Eq. (jS)) 
and expressing k^a^ = K/Ef, and k^a'^ = K/E'^, we 
derive a simple and very useful form of the lineshape, 
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In the limit fc — > 0, where 5' = —ka', we have 
Ff ~ K^/"^, which is the Wi gner threshold regime. The 
transition rate peaks between if ~ for |a'| ^ a and 
K = Eh/S for \a'\ <C a and decreases to higher K. For 
very large K > It? /mr^, Eq. © is no longer valid. A few 
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The unit Franck-Condon factor for a' < means that 
within the approximations we made, the molecular wave- 
function (pmir) is fully expanded by the scattering states 
\K). For a' > 0, we have Sf < 1, which implies there 
are additional outgoing channels. This is because the 
additional bound-bound transition process is opened for 
a' > 0. 



B. Bound-bound transition 

For positive scattering length in the outgoing channel 
a' S> > 0, a weakly-bound molecular state |m) ex- 
ists. Bound-bound transitions are allowed. From energy 
conservation, we have 
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Effectively, K = — i?b- The bound-bound Franck-Condon 
factor Fh{K) can be calculated similar to Eqs. ^ and 
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where we have used the molecular wavefunction in Eq. © 
for both |m) and \m') states and have introduced the S- 
function to provide energy normalization analogous to 
that of Ff{K). The bound-bound line strength is then 
Sb = 4aa'/{a + a')^. 

Note that sum of the bound-bound and bound-free in- 
tegrated line strengths is identically one for a' > 0, as 
expected from the wavefunction projection theorem: 



S*/ + S'h = 1 . 



(16) 



III. 



(10) 



COMPARISON WITH NUMERICAL 
CALCULATION 



We check the validity of the above analytic formulas by 
comparing them to a numerical calculation. We choose 
fermionic ^Li atoms as our model system, for which the 
interaction parameters are precisely known by fitting var- 
ious cold atom and molecule measurements to a multi- 
channel quantum scattering calculation Q. The energy 
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FIG. 1: Energy structure and the scattering lengths of ^Li2 
in the (1,3) and (1,2) channels (dotted lines). All energies 
are referenced to the (1,2) scattering threshold. In these 
two channels, Feshbach couplings induce the formation of 
molecules (solid lines) below 69.1mT and 83.5mT, respec- 
tively. Arrows show the bound-bound and bound-free tran- 
sitions based on molecules in the (1,2) channel. The lower 
figure shows the scattering lengths in the two channels. 



structure and the scattering lengths in the two relevant 
channels A = (1, 2) and A' — (1, 3) are shown in Fig. ^ 
Here, (1, 2) refers to the state with one atom in state |1) 
and one in |2), and \N) refers to the lowest A^th internal 
state in the ^Li atom ground state manifold. 

To calculate the scattering phase shifts, we construct 
a reduced single-channel Hamiltonian to describe the 
(1,3) continuum. The phase shifts from this model are 
nearly indistinguishable from those from the full multi- 
channel calculations in this range of B-fields. The van 
dcr Waals interaction is set to Cg = 1393.39(16)a.u. 
0, yielding an interaction range of — 29.884(3)ao. 
The scattering lengths are set to match the values from 
the multi-channel calculation. The resulting scattering 
phase shifts from the calculation are compared with the 
effective range expansion given in Eq. © (see Fig. 
Here, the scattering parameters are based on the (1,3) 
scattering states at 72.0mT with a scattering length of 
a' = — 7866ao. The result shows the expected behav- 
ior that fccot(5' = — 1/a' provides an excellent fit at low 
scattering energy i? < fc^ x 20/i K and the effective range 
correction works up to < fc^ x lOmK. Here fc^ is Boltz- 
man's constant. 

We also compare the Franck-Condon factors obtained 
from the numerical calculation and from Eq. H1U|) . The 
numerical calculation is based on an initial bound state 
in the (1, 2) channel and the final state (1, 3) continuum, 
described by the reduced Hamiltonian. We show the cal- 
culations at two magnetic fields of 72.0mT and 68.0mT, 
where the scattering parameters are set according to the 
multichannel calculations. The results show that Eq. H1(J|) 
provides an excellent fit to the full bound-free spectra, see 
Fig. El 
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FIG. 2: Comparison of the scattering phase shifts from nu- 
merical calculation and from the effective range expansion. 
The numerical calculation (open circles) is shown together 
with the predictions from Eq. © with (solid line) and with- 
out (dashed line) the effective range correction. The scatter- 
ing parameters are based on ®Li in the (1,3) channel near 
72.0mT with aw = -7866ao. 
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FIG. 3: Comparison of the Franck-Condon factors from nu- 
merical calculation and from theory. Numerical calculation at 
72.0mT (open squares) and 68.0mT (open circles) are plotted 
together with the formula Eq. HIOII (solid line at 72.0mT and 
dotted line at 68.0mT). The parameters are Eb = ksX 16.2^K 
and a' = 25173ao at 68.0mT; Et = x 6.2^iK and a' = 
-7866ao at 72.0mT. 



IV. RF LINESHAPE NEAR A FESHBACH 
RESONANCE 

The appearance of the bound-bound transition for only 
a' > seems to suggest a distinct behavior in rf excitation 
near the Feshbach resonance where the scattering length 
changes sign. In this section, we show that the evolution 
of the rf spectrum is actually continuous when a i?-field is 
tuned through the resonance position. We use Eq. H1U|) 
andEq. H15|) to show the bound-free and bound-bound 
spectra in the vicinity of a Feshbach resonance. 

Figure ^ shows a continuous change near a Feshbach 
resonance from a bound-free lineshapc for a' < to a 
combination of bound-bound and bound-free transitions 
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FIG. 4: Transition from a bound-free transition to a bound- 
bound transition with Ei, = 1. For the bound-free transition, 
the Franck-Condon factor is plotted as a function of K/Ei,. 
Curves of different a' /a are offset and shown in the order of 
that near a Feshbach resonance in the A' channel, (dotted 
line for the resonance condition) For a > 0, the locations 
of the bound-bound transitions are indicated by the sharp 
peaks. Notably, bound-free transition vanishes at a' = a. 



for a' > 0. When the scattering length approaches neg- 
ative infinity, the linewidth of the bound-free lineshape 

^ and "evolves" into the 
^' > 0, 
a, the 



approaches zero as ^ E"^, ^ a' 
bound-bound delta function. Remarkably, for a' 
the bound-free transition is much weaker. At a' - 
bound-free transition is fully suppressed, ~ 0. 



The vanishing bound-free transition for a' ~ a can be 
understood by the wavefunction overlap given in Eq. |0J. 
As a positive scattering length a' indicate a zero in wave- 
function near r « a, the sign change of the scattering 
wavefunction results in a cancellation in the overlap. Al- 
ternatively, we may consider the transition amplitude to 
the continuum A' is exactly cancelled by the Feshbach 
coupling between molecular state |m') and the contin- 
uum A' . This interference effect results in the Fano-like 
profile of the peak transition rate near a Feshbach reso- 
nance 01 (see Fig. El • 

On the other hand, the bound-bound Franck-Condon 
factor approaches one when a' « a, since both molecular 
wavefunction are identical for r > tq. The peak bound- 
free transition rates therefore shows a strong asymmetry 
with respect to the sign of the scattering length, shown 
in Fig. 

Based on the above features, the rf spectroscopy does 
provides a new strategy to extract various important 
scattering parameters. In particular, while cold collision 
measurements are generally insensitive to the sign of the 
scattering length, our rf excitation spectroscopy is dras- 
tically different for the two cases and provides a new tool 
to probe the scattering properties. 
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FIG. 5: Fano profile near a Feshbach resonance in the outgo- 
ing channel A' . The peak Franck-Condon factors with Ei, = 1 
are calculated for different a /a based on Eq. HUH . Near the 
Feshbach resonance, a/ a' approaches zero. 



V. COMPARISON WITH EXPERIMENT 

To compare with the experimental result, we take ^Li 
as a model system and calculate the rf spectra. The 
experimental spectra are obtained from the Li group in 
Innsbruck 

Adopting the convention used in Ref. |2| , we define the 
rf offset energy as _B = Enp — Eq = K + E^. Using 
Eq. H1(J|) . we can rewrite the Franck-Condon factors as 
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-5(E -Eb + E'b) for a' > (18) 



We compare the theoretical curves with the experi- 
ment in Fig. 13 An integration bandwidth of IkHz is 
assumed. When the magnetic field approaches 69.1mT 
from high values, the molecular binding energy Eb in the 
(1,2) channel increases (see Fig. This dependence 
is shown in Fig. |B1 as the whole excitation line moves 
toward higher frequency. At the same time, the line- 
shape becomes sharper due to the Feshbach resonance in 
(1, 3) channel at 69.1mT. Below 69.1mT, we expect both 
bound-free and bound-bound transitions are allowed. At 
68.4mT, these two components coexist, but cannot be 
clearly distinguished due to the small binding energy of 
E'fj = h X 0.5kIIz, lower than the experimental resolu- 
tion of IkHz. At 67.6mT and 66.1mT, the bound-free 
transitions are strongly suppressed. The bound-bound 
transition then shows up as the dominant component in 
the spectra. The excellent agreement between the exper- 
iment and the calculation over a large range of magnetic 
fields is remarkable. This result strongly suggests that 
the rf molecular spectroscopy can be used for a precision 
determination of the atomic interaction parameters 
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FIG. 6: RF spectra of ^Li2 molecule at different magnetic 
field values. Theoretical curves are based on Eq. I17II and 
Eq. If 81 with the binding energies and scattering lengths ob- 
tained from the multichannel calculation. To shown both the 
bound-free and bound-bound spectra, an integration band- 
width of IkHz is assumed. RF-induced loss on the molecular 
population at 76.4mT (solid squares), 72.0mT (open squares), 
69.5mT (solid circles), 68.4mT (open circles), 67.6mT (solid 
diamonds) and 66.fmT (open diamonds) are re-scaled. Ex- 
periment data are from the Grimm's group in Innsbruck. 



VI. CONCLUSION 



the recent rf measurements on Li2 molecules. 

An interesting case is studied when a Feshbach res- 
onance occurs in the outgoing channel. We show that 
the bound-free spectra in the absence of a bound state 
smoothly evolves into a combination of bound-bound and 
bound-free spectra when the bound state formed near 
a Feshbach resonance. The bound-free transition rate 
strongly depends on the sign of the scattering length and 
shows a Fano-like structure near the resonance. 

We like to point out that the rf spectroscopy based on 
weakly-bound molecules can be an excellent tool to de- 
termine the long-range interaction properties with high 
precision. From the excitation spectra, the molecular 
binding energies, the atomic scattering lengths and their 
signs, and the scattering phase shifts can be determined. 
Furthermore, in contrast to conventional cold collision 
measurements, rf transitions on molecules are insensitive 
to the sample density and temperature, can probe differ- 
ent scattering channels by tuning the rf frequency to dif- 
ferent states, and can be implemented instantly without 
thermalization. Finally, rf transitions between molecu- 
lar states may also provide a new avenue to transfer the 
molecular population to low-lying molecular states. 
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